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I . INTRODUCTION 

Recent interest in STOL aircraft has stimulated investigations of 
airfoil configurations employing both powered and unpowered devices for 
obtaining high lift coefficients and improved stability characteristics 
at low flight speeds. The practical feasibility of the powered schemes, 
such as jet flaps, blown flaps, and augmentor wings was greatly enhanced 
with the advent of high by-pass ratio engines. This provided substan- 
tial incentive for obtaining theories capable of predicting induced 
velocity fields and forces for configurations employing such devices. 

A jet flap comprises a thin jet of air exhausting from the trail- 
ing edge of the airfoil at an angle to the direction of the freestream 
(refer to Figure 1). The presence of the jet contributes to the lift of 
the airfoil in two ways: (1) the aerodynamic field is altered by the 

jet, increasing the net lifting pressure force on the wing’s surface, 
and (2) the reaction to the momentum flux of the jet contributes to the 
lift. Two related devices, the blo^sn flap and the augmentor wing are 
also shorn in Figure 1. In the blom flap scheme the high velocity jet 
is ejected over the upper surface of a mechanical flap utilizing the 
Coanda effect to direct the jet streaju. Tlie augmentor wing employs an 
ejector system combined \d.th the trailing edge flap design in order to 
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aqgpaent the thrust of the primary jet. This latter device develops an 
increased thrust and also produces smaller nose-dcmn pitching moments 
than the other devices. The jet flap remains the simplest of the 
configurations . 

One of the first analytic investigations of the jet flap schane 
was a two-dimensional theory developed by Spence.^ This was later ex- 
tended to handle the three-dimensional problem with lijiuLtation to 

elliptic loadings, constant jet momentum coefficients, and constant jet 
• 2 

deflection angles. Several theories seeking to provide more generality 

have been published. Those of Kemey and Tokuda^ who employed matched 

asymptotic expansions to provide a method capable of handling more general 

wing-flap configurations and the lifting surface theory of Lopez and Shen^ 

which offered the further advantage of also providuig solutions for 

augmentor v^dngs and blow flaps are representative. Other theories which 

have considered the ground effects problem^ or provided more rapid compu- 

tation techniques are now available. All of these methods, however, 

employed approximations inherent in linear theories. That is, all flow 

angles were assumed to be small. Although it might be expected that 

such theories should be capable of providing good results for the forces 

and moments on the airfoil, at least for situations where the jet flap 

is at a small angle to the free stream, or where the momentum flux is 

small, the accuracy with which the downwash field could be obtained is 

open to serious question. In the past few years, two-dimensional 

8 9 

theories wliich accounted for nonlinearities have been reported, * 
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Each o£ the existing theories may adequately account for parts o£ 
the overall problem, but none may be expected to be capable o£ treating 
the entire spectrum o£ flight regimes and jet loadings. With low jet 
loadings and sufficient airspeed, for example, the field induced by (the 
jet may be regarded properly as a small perturbation to the basic' aero- 
dynamic field. With higher jet loadings and lower airspeeds, this may 
not be acceptable. If computation of the downwash field near airframe 
manbers such as the tail surface is to be required, along with forces and 
moments on the wing, additional care in the treatment of the jet is re- 
quired. And there are regimes of flight demanding yet other factors, 
such as consideration of the ground plane, be included in the theory. 

The analysis described herein is intended to liandle the situation 
v^ere the induced aerodynamic field of the jet is not small in comparison 
with the normal aerodynamic field of the vehicle. It is to apply to 
wings with jet flaps of either high or low loadings and offers the 
additional benefit of improved downwash computations. 

Formulation of the complete inviscid, incompressible, three-dimen- 
sional jet flap problem is formidable. Even if the thin sheet approxi- 
mations are employed, it is still necessary to apply a nonlinear boundary 
condition at the jet boundary, the position c£ which must be determined 
in the solution. It is essential, therefore, to employ further simplifi- 
cations in order to obtain a solution for the nonlinear three-dimensional 
problem. It was considered important to account for the position of the 
vorticity associated with the jet, departing substantially from tlie plane 
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of the airfoil in the flight regime of interest. Mien rollup of the 
jet sheet may be neglected, as it is here, it is necessary to introduce 
either an approximation for the spanwiso geometry of the jet or the 
functional form of the spanwise distribution of vorticity in the jot, 

Ihe latter offers significant simplification of tlie theory and was chosen 
in this theoretical model. The streamwiso variation of the deflection 
of the jet sheet remains to be solved as an iiiqiortant element in the 
solution. Ihe specific assumptions and approximations made in tlie pre- 
sent theory, are discussed in Chapter II along with a general discussion 
of the solution procedure. Cliapters III, IV, and V are concerned with 
the detailed development of the theory. The theory is then applied to 
a number of situations and compared with a few of the linear theories as 
well as the experimental results of Williams and Alexander^^ in Chapter 
VI. Chapter VII concludes with a discussion of the strengths and weak- 
nesses of the theory as well as areas of possible further study. 

II. FORMULATION OF THE ANALYTICAL MODEL 

A theoretical model attempting to represent the forces, moments, 
and induced aerodynamic field in the vicinity of an airfoil equipped 
with a jet flap must represent the jet reasonably well within a few 
chord lengths of the wing. This was an accepted premise of the analysis. 
^ virtue of the dominant streamwise momentum flux of the jet in the 
near field, the rolling up process of the streamwise vorticity in the jet 
sheet plausibly could be neglected in the near field region. But the 
position of the jet sheet was retained as an important element in the 
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theory even though the precise spanwise geometry would not be represented. 
It is this feature which primarily distinguishes the present theory fran 
the linear theories. Viscous effects are ktiom to contribute to the 
induced field of the jet, but they are neglected along with compressibility 
effects which play a role in the flow ivithin the jet. The entrainment 
properties of the jet, however, could be added to the theoiy with little 
difficulty. 

Even neglecting viscous and compressibility effects from the outset 
in contemplating a theoretical model for a jet flap, an intractable, three- 
dimensional, inviscid problem in aerodynamic theory remains. It tvould be 
necessary to obtain solutions for the flow within the jet (the internal 
flow problem) and the external aerodynamic field, subject to the pressure 
(dynan\ic) and kinematic boundary conditions at the interfaces between the 
two flows. Tliese interfaces could be expected to be highly convoluted 
sheets of semi-infinite extent. The geometry of the sheet is not known 
initially and iviould have to be obtained as part of the solution. It is 
essential to introduce further assumptions to obtain a tractable problem 
with an appropriate reduction as to what may be expected of the theory. 
Moreover, for the information sought in the analysis (forces and moments 
on the airfoil and downwash data), a precise solution of the complete in- 
viscid problem would seem imnecessarily detailed and unwarranted. 
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A. Principal Assumptions 

The principal assumptions in the formulation of the analytical 
model may be outlined c'onveniently as follows. 

1. The flow field is assumed to be inviscid and incompressible. 

The inviscid approximation, as in most aerodynamic theories, sacrifices 
tlie capability of accounting for boundar>' layer phenomena such as skin 
friction on the airfoil surface and anticipation of separation phenomena. 
These phenomena may be treated in a separate manner, given the basic 
aerodynamic flow field, as in more conventional aerodynaunic problems. 

IVfo aspects of the inviscid approximation relevant to the jet flow which 
could be of some importance in the jet flap problem are the viscous en- 
trainment by the jet and the growth of the jet thickness by virtue of 
mass and momentum exchange with the external flow. Ihe former is likely 
to be most importajit. The aerodynamic field induced by viscous entrain- 
ment could be accommodated by employing a suitable distribution of sources 
and sinks in the inviscid model. This is left for possible extension of 
the present theory. Under the conditions of higlily loaded flaps, the con- 
tribution should be a sinall one. Compressibility may be of significance 
in some pdrts of the jet flow and possibly in the concentrated vortex 
cores of the wake. Considering the fact that most of the jet is in a low 
subsonic flow regime and that the near field is of more importance than 
the distant wake, compressibility should not exert a significant influence 
on the results. For higher freestream velocities where the aerodynamic 
field might be locally compressible, compressibility could possibly be 
accounted for by employing the Goethert rule. 
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2. The airfoil and the jet flap are assumed to be symmetric about 
a plane containing the undisturbed frees tream velocity vector. This is 
an essential part of the analysis; more general geometries would recpaire 
a much more elaborate analysis . 

3. Airfoils equipped with both full span and part span jet flaps 
with constant jet momentum coefficients and which issue from a straight 
trailing edge are considered. Variable distributions of the jet momentnnt- 
coefficient would be admissible in the basic theory, but such situations 
were not considered in the computations. 

4. The thickness of the wing is neglected; however, camber and 
variations in the wing planform are included. This procedure follows the 
usual -thin airfoil approximation in aerod>Tiamic theory. The resulting 
force system on the airfoil may later be mnriified to account for the 
thickness effects of the airfoil in an approximate manner by mploying 

a modified lift slope as initially used by Spence^ and discussed in 
greater detail by Lissaman. Furthermore, such effects also miglit be 
represented by a suitable distribution of sources along the surfaces 
used to denote the airfoil. 

5. Several approximations are introduced pertaining to the repre- 
sentation of the jet. 

Ca) The thickness of the jet is neglected, and the jet momen- 
tum coefficient is taken to be constant. Interaction between the jet 
and the external field then arises solely from the normal force due to 
pressure differences across the jet sheet. The jet is assumed to have 
no spanwise con^onent of velocity. 
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(b) Rollup of the jet sheet at its edges is not considered. 

TMs may be of importance for lightly loaded jets, at high freestream 
velocities, or far enough downstream. For the situations of interest 
here, however, the rollup is likely to be delayed to a point far enough 
downstream to be of little influence on the results. 

Cc) The vorticity associated witli the jet is assumed to lie in 
a sheet generated by a system of straight lines normal to the free stream, 
direction (refer to Figure 2a) and tangent to the actual jet sheet at the 
plane of symcnetry. The position of the .stream surface associated with the 
jet may not be composed of straight line generators, as indicated in Figure 
2a, but this arrangement does place the vorticity close to the actual 
surface, an important factor Avhen the jet lies significantly below the 
plane of the wing. Ihe small difference between the actual jet sheet and 
the sheet on whicli the singularities are placed should not comproiaisc 
cue solution in a serious manner. Similar approximations have mot with 
ample success in conventional aerodynamics'. 

(d) Tlie lateral (spanwise) distribution of the sp<mwi.so com- 
ponent of vorticity in the jet sheet is assumed to be elliptical. Ihis 
distribution falls to zero at the edges of the -jet in the proper manner 
and greatly simplifies the theory. To be sure, the spanwise distribution 
in an actual jet will not follow this pattem, but will vary in the 
.streamwise direction leading to concentrated vortex cores far doMistream. 
Nonetheless, in the near field of primary importance, at least, the 
distribution probably will have a form approximating an elliptical loading. 
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Moreover, in view of the neglect of the fact that rollup occurs, a plausible 
assun5)tion of this sort remains consistent with the framavork of the 
balance of the model. It should be observed that employing such an as- 
sumption for the spanivise distribution of the vortex strength is tanta- 
mount to specifying the lateral geometry of the jet and it is far more 
effective in simplifying the theory . One or the other must be given if 
rollup is neglected. In order to provide greater care in modeling the 
jet vMle retaining convenient working equations, the former was chosen. 

Hie merit of this assumption ultimately may be tested, of course, by 
comparisons with experimental data, 

6. Certain assumptions are made as to the mathematical conver- 
gence of the solution procedure. This method amounts to (a) starting with 
a first approximation for the flowfield of the wing alone from an ele- 
mentary linear theory, (b) solving the jet problem subject to the free- 
stream and the disturbance field of the wing alone, and (c) employing 

the jet solution so obtained, together with the undisturbed freestream, 
as the externally Imposed field in a more complete solution for the 
disturbance field of the wing. If these last results prove unsatisfac- 
tory, the procedure may be continued, obtaining a new jet flow, etc. 

It is assumed this procedure does converge to a solution. 

7. In computations involving the influence of the disturbance 
field of the airfoil, the vorticity in the wake of the wing alone is 
assumed to lie in the plane of the surface used to represent the airfoil 
in a manner sinular to existing lifting surface theories. For a lightly 
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. loaded wing and in the region close to the airfoil, this vorticity \>dll 
have a primary effect of distorting the jet sheet laterally. Because 
this effect is not a primary consideration, and since the influence of the 
disturbance field of the wing is sm^l compared with that of the free— 
stream for the jet solution, such a contribution night be reasonably- 
neglected altogether. However, the wake vorticity is retained and treated 
in the above manner because it does restilt in simpler equations. A more 
accurate accounting of the vorticity shed by the wing would add signifi- 
cant con^lication to the analysis unless some assun 5 )tions as to its 
position couLLd be made a priori. That part within the spanwise extent, 
of the jet presumably would follow the jet con torn:. But that near the- 
wing tips would lie somewhere between the plane of the vring and the 
location of the jet for an airfoil equipped with a part span jet flap. 

B. Mathematical Foimulationy 

The coordinate systems employed in the analysis, along with tha- 
principal geometric parameters of concern are illustrated in Figure 2. A 
right-handed cartesian coordinate system is employed, taking the x-axis 
alcmg the direction of the freestream velocity, the y-axis alcnig tlie 
intersection of the airfoil's trailing edge and the jet, and the z-axis 
in the vertically upward direction. For some confutations in the jet 
analysis it is convenient to employ an orthogonal curvilinear coordinate 
system located id. thin the jet sheet: ^ in the streairodse direction, n 
the spanwise direction, and c normal to the sheet. 
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Referring to Figure 2.b, the two principal angles are defined in 

the cross sectional view. The angle of attack C«) is the angle measured 
from the x-axis to the chord line and the jet deflection angle (t) from 

the chord line to a line which is tangent to the jet as it leaves the 

trailing edge of the airfoil. The notation pertinent to the span lengths 

and the wing chord are also shown in the figure.. 

The positions of the physical elements employed in the problem, 

the jet sheet and the (thin) wing, are sho^vn in bold outline in Figure- 

2a. The locations of the singularities used to represent these elements 

are shorn in lighter outline. The singularities used to represent the 

airfoil surface are placed in -the z = 0 plane in an area corresponding 

to the planform of the wing. Those for the jet are positioned on a 

surface tangent to the actual jet sheet in the plane of symmetry (the- 

r 

y * 0 plane) having straight line generators (parallel to the y-axis) 
as already mentioned- 

Given a specific \\ring planform and camber distribution, together 
with the span length of the jet flap, the necessary parameters for this:v 
problem aret 

(a) the angle of incidence of the wing (a) 

(b) the frees tream velocity (U^) 

(c) the angle of deflection of the jet (x) , and. 

(d) the jet momentum coefficient C^j). 

The objective of the analysis, then, is to obtain results for the forces 
and moments on the airfoil, and the velocity field at points of interest 
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with respect to the airframe design (e.g. , the downwash in the vicinity 
of tail surfaces). 

1. General Relations 

Witliin the framework of the assumptions discussed above, the 

analytical problem amounts to one of the Neumann type in’ potential theory. 
The flow is everywhere irrotational and solenoidal* giving rise to a 

representation of the field in terms of a scalar potential. The normal 

derivative of the scalar potential on the airfoil and on the jet sheet 

must be zero; these are stream surfaces. The great difficulty of the 

problem arises from the fact that the position of the jet sheet is not 

known initially and must be determined in the solution. This requires 

matching the pressure difference across the sheet with the normal 

acceleration of an element of fluid in the j et (dynamic condition) . 

The relations required in the theory may be obtained from the 

ll 12 

Fundamental Theorem of Vector Analysis (Helmholtz's Theorem) ’ . Re- 

ferring to Figure 3a, consider for the moment a closed, singly connected 
domain surrounding the wing and jet and containing both (the domain en- 
visioned here is large enough to contain even the starting vortex of the 
system) . Define the perturbation velocity v in the usual maimer 

V = V - e^ 

where ^ is the complete velocity vector. Then 

V = -V(J>.+ V X X 

Excluding the singular surfaces representing the wing and the jet. 
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where 


and 


■ w III 
III 


71 . 
R 


dT' 


V’ X v‘ 
R 


dT’ 


R = 1^1 = y(x-x’? + (y-y*)^ + (z-z’)^ 


In these relations, integration is with respect to the source coordinates 

(x') and the potentials (j) and ^ are functions of the field point x. 

Combining the above equations a general relation may be obtained for the- 
■^13 

vector V 



Ao /S> y 

S represents the stxrface bounding the region T and n represents the unit 


outward normal of S* 

The vector field of concern is solenoidal (V 
therefore, the first term of Equation (1) is zero, 
enough, v is zero at every point on the surface S. 
is entirely given within the domain as 


• V = 0 everywhere) r 

A, 

For a domain T large 
Then the vector fieldi 
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This situation, of course, is not quite what is needed for the 
problem at hand. Rather, one is interested in a part of the domain in 
which the floi^ and jet geometry are stationary in time, a domain con- 
tained \Nd thin the one described above and far enough removed from the 
time varying starting vortex region so as not to be influenced by it 
(refer again to Figure 3a) . From the foregoing discussion, it should be 
clear that the velocity at a point of interest within the inner domain 
will be determined as the sum of two integrals of the form of Equation (2) 
one over the inner domain and one over its complement. Since these re- 
gions may be chosen for our purposes to be as large as desired, the con- 
tribution of the second integral to field points of interest within the 
inner domain (not too near the intersection of the jet sheet with the 
boundary of the domain) can be made small enough to be inconsequential 
in the analysis. Equation (2) , then, may be written as 


v(x) = 



(V* X V') X ^ 

IT3 


Now, the vorticity in the domain of interest is assumed to be con— 
centrated in the form of vortex sheets. These singular sheets are re- 
lated to the actual distributed vorticity regions by introducing the 
vortex strength of the sheet, y» such that (refer to Figure 3b) 

lim m t dS = Y dS 
t 0 
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H^iere w is the vorticity vector (to = V x v) ajid the element of area S 
lies in a vortex region of thickness t. It may be shom^"^*^^ that the 
strength of the vortex sheet is related to the vector difference of 
velocities on either side of the sheet. 

Y X n = A V = (v^ - v~) (3) 

where n is a unit noi-mal to the slieet and v (v. ) is the velocity v'^ector 
above (below) the sheet. Employing the definition of y» the velocity 
field v(x) may now be expressed in the form 

where S now represents the totality of vortex sheets within the domain^ 
of interest. Equation (4) is the starting point of this analysis as it 
is for most other problems in aerodynamics- 

The problem now is to determine the vortex strength distribution , 
on these sheets so that the flow is everywhere tangent to them. Suppose 
for the mcanent the sheet positions may be described in the fomt- 

f (X, y) - z » 0 

The kinematic condition of tangency then becomes: 

^ • 7 (f - z) « 0 


or 
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Furthermore, since y must be solenoidal and since by definition it lies in 
the vortex sheet, it has only one independent component. Equations (4) 
and (5) then provide four equations for the variables u, v, w and, say,. 

Yy. If the positions of the vortex sheets were known, then, by substitur- 
tion of Equation (4) for u, v, and w, along with a relation between y^ 
and Yy, obtained from the solenoidal character of Yi into Equation (S) , a 
single integral equation over a two-dimensional sheet results for the 
determination of Yy* In view, of the assumption to be made as to the 
spanwise distribution of vorticity, the condition given by Equation (5) 
mav be met at the plane of symmetry where the sheet of singularities 
and the stream surface of the jet coincide. Application of Equation (5), 
starting at the trailing edge of the airfoil, may then yield the form 
of the actual jet sheet at points not in the plane of symmetry. This will 
not, in general, coincide with the sheet of singularities representing 
the vorticity of the j et (refer to Figure 2) , 

The location of the vortex sheet associated with the jet un- 
fortunately is not kno\\n ini.tially. ^fc)reover, the jet sheet must be 
positioned such that the pressure force acting on it balances the local 
normal acceleration of the flow within the jet. This dynamical condition, 
then, provides the necessary additional relation for the determination of 
the position of the jet sheet.- 

Consider a jet sheet generated by straight lines such that it has 


only one finite curvature and with its normal vectors lying in planes 


parallel to the x-z plane. i\pplying tlie momentum theorem o£ fluid 
mechanics to an element of the jet sheet, it may be shown simply that 

+- 


K 



where k represents the curvature of the jet sheet, p^ Cp") is the pres- 
sure above (below) the sheet, and J is the momentum flux of the jet per 
unit length of span. This relation may be placed in a more suitable form 
as follows. By employing Bernoulli's equation along with Equation (3) 
the pressure difference across the jet may be written 


AP-p"*"-? - - j p - V"^) * -p (y X n) • Vjjj 

where 


is the mean velocity at a point located on the jet sheet. Utilizing the- 

2 . 

definition of the jet momentum coefficient 


r 

I 


_j 


the dynamic condition may be written in the form 

_ 2 Cy X n) • 

K ® 2 

Cj u: c 


( 6 ) 
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Bjr definition, the curvature is related to the geometry of the jet sheet 
under consideration through tlie relation 



d% 

S? 



dx 

|2j3/Z 


vihere the surface in question is given by the relation 


f(x) - z » 0 

This form of the dynamical condition is to be satisfied at points on the 
jet sheet in the plane of symmetry where tlie jet sheet coincides with 
the sheet of singularities representing the jet vorticity. This condi- 
tion is not met strictly at any other point on the jet sheet, of course, 
within the approximations involved in the analysis. This is again a 
consequence of the neglect of rollup in the theory. 

liquations (.4J, (^5J, and (.6j provide the basic relations from which, 
the vortex strength distribution and the position of the jet sheet may be 
detezmined. An iterative procedure was adopted for that purpose in the 
analysis. 

2. Method of Solution 

The solution procedure adopted for the analysis is most conveniently 
presented by considering the perturbation velocity field on the airfoil 
and the jet separately. Ert5>loying equation (4) the required relations 
■ay be e:q>ressed as 
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( 7 ) 


7cxj) - 
V(x^) =• 



IT 

xul ds- t 

T.-3 



S-ilds' 


Y' X ^ 

r3 


dS» 


(8) 


where the subscript J (A) refers to the jet (airfoil). Equations (7) and 
(8) provide a set of coupled equations for the vortex strengths on the 
airfoil and the jet. They are analogous in some respect to the bi-plane 
problem. In choosing a method of solution for this problem, it is helpfi.il 


16 


to keep ill mind the physical conditions of primar>’ interest. Tlic influence 


of the disturbances produced by the wing on the jet (represented by tlie 
first surface integral in Equation (7)) may be expected to be small com- 
pared with the interaction of the jet and the freestream for the flight 
conditions of interest, idiereas the reverse is not true. The jet induced 
disturbances (represented by the second surface integral of Equation (8)) 
should have a significant influence on the wing field for a highly loaded 
jet. 

A solution for this set of equations may be obtained as follows. 

(a) A solution for the wing alone is obtained as a first approxima- 
tion for the use of modeling the influence of the wing on the jef in 
Equation (7) . This may be as simple as a lifting line approximation or 

a more detailed solution if readily available. 

(b) A solution for the jet is obtained subject to the influence of 
the combined fields of the wing and the freestream. This solution com- 


28 


I 



prises obtaining the distribution of the vortex strength along the jet 
and the position of the jet sheet 

(q) The j et solution and the freestream are taken as a new in^josed 
field for a detailed solution of the wing problem. The solution pro- 
cedure employed herein is of the lifting surface type^.. 

Given these results, the forces and moments on the wing may be- 
obtained by direct integration, and the velocity field at points of 
interest in the vicinity of the wing or jet may be computed as needed. 
Should the results for the last wing solution prove ijiaccurate, it is 
possible to employ that solution as a new initial approximation and 
repeat the entire procedure. 

The analysis of the jet subject to the influence of tJie freesti'onm 
and the disturbance field of the airfoil mdy be referred to here as the 
jet problem . ITiat of determining the aerodynaiaic field of the wing sub- 
ject to the freestream and the field of the jet may be termed the airfoil 
problem in the following treatment. 

3. Working Equations 

A solution for the jet problem requires deteimination of the dis- 
tribution of vortex strength (Yj)-on the jet sheet as well as the loca- 
tion of the jet in the plane of symmetry (z = fj(x,0)). For a given fj, 
obtained either by £in initial assumption or from the results of a previous 
iterate, the vortex strength may be determined by satisfying the kine- 
matic condition, 

[Uco \ vCx,0,fj)] . V Cfj - z) = 0 (9) 
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\ 4 here 




( 10 ) 


and accounts for the influence of the do\vnwash field created by the 
airfoil on the jet. The proper solution is one \vhich also satisfies 
the dynamic condition 


K 


^ ^ n,) • i 


(U) 


An iterative procedure is adopted for varying £j(x) so that Equation (11) 
is met in the solution, along with Equations (9) and (10), 

In a similar manner for the airfoil problem, a solution for the 
vortex strength distribution of the airfoil (y^^) may be fomid by applying 
the kinematic condition 


C^co ®x ^ ■ 2) “ 0 


( 12 ) 


where 




to the location of the projection of the airfoil on the z 
the above equations 


(13) 

0 plane. In 


-2 = 0 
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represents the surface of the airfoil and Vj accounts for the dowro^ash 
field created by the jet on the airfoil. 

III. JET ANALYSIS 


Consider now the problem of determining the vortex strength distri- 
bution and the location of the jet sheet, subject to the influence o£ 
both the freestream and the disturbance field created by the airfoil 
For present considerations, the field induced by the airfoil will be 
assumed to be known. 

Recalling the curvilinear coordinate system (5,n) located in the 
jet sheet, the vortex strength distribution associated with the jet may' 
be expressed as 


YjK.n) - Yj 5 % * Yj„ % 

Since the vorticity field is solenoidal, and Yj^ are related through, 
the relation 


55 5n ■ 


or 


'Jn 


® Yjr 


-s, 


The assumption that the spanwise distribution of loading on the jet has 
the same functional form at all streamwise positions along the jet implies 
that the con^onents of the vortex strength may be represented in terms 
of products of the form 
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then 


= G(5) HCn) 


Yjn “ - ^ 


Defining 


h(Ti) = H(n') dn’ 


the equation for Yj^ be written 


3 GCg) 


If the loading on the jet is to be elliptic, as discussed earlier. 


and then 


h(n) = 


H(n) = 


1 - Cn/sj)' 


1 - (n/sj)‘ 


The choice of an elliptic loading prov^ides the primaiy advantage of 
simplifying some of the integrals in the subsequent elements of tlie 
theory, and does have the proper characteristics at the etlges of the jet. 
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A. Integro-Differential Equations for the Jet Problem 


1. Analytical Relations 

Recall first the relaticxis given by Equations (9) to (11) 


( 16 ) 


(Uott + u) - W - 0 


C17) 


K = 


d^f. 


dx 


2 CYj X nj) . % 


m 


Fra 


j\2]3/2 


Cj V c 


(18) 


where these relations all pertain to field points located on the jet sheet 
in the plane of symmetry (y = 0) . Then 


y = 0 


= (U„ + u) + w e^ 


since v =' 0 in this plane . And 




It has been determined that integrations carried out with respect 
to X and y rather than £ and n are most convenient. Ihe cartesian para- 
meters involved may be expressed as 
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TjX “ 'Iji » ' 


''Jr ’ ’fJn 


w(‘ • (Sir 


df, 


yjz “ yjK ® - a^/[i ♦ 1 ^)^''^ 

dSj' = d?’ dn' = sec 0' dx' dy' 

ffj ’ * \ ■ -^ j / - -^Jz' O ' - y ) 


(19) 

( 20 ) 
( 21 ) 


where 


CYj ’ X R ) j , » Yjx ’ (y - y ') - Yjy'(x - x») 


R = {[X - x']2 - [y - y']2 + [z - f, (x*)]^} 


1/2 


and, referring to Figure 3 

sec 0 


vMSF 


Then the components of Equation (16) take the form 


u(x.0,z) - u^(x,0,z) 


sec 0* dx* dy' (22) 


w(x,0,z) = w^(x,0,z) + ^ 


7 -Yjx' y' “ Ttv' 


“T 


'Jy 


- X’) 


sec 0* dx' dy' (23) 


-SjO 
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Introdi’cing Equations (14) and (15) for Yj^ and Yj^ into Equations 


(19) through (21), 


Y- = G H/ 1 + 


/r /df \2il/2 

H/hMJ 


and let 


-ar 


Y-t, - G 




g(x) = G [5(x,fj)] 


TTien since ^ ^ ^ ^ cos 0 , Equations (22) and (23) may be written 


u(x,0,z) = u^(x,0,z) 


!_ I I >5^ 


1 "j' ? 

g' T^dx' 
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~ J 

w(x,0,z) = w^(x,0,z) - ^ I g' dx' j 

o -St ^ 


W / ^ 


(x - X') dx' 


J 

f h' dy' 

“F" 


( 25 ) 


The integrals with respect to y' may be evaluated directly (refer to 
y^pendix A) to yield 

df * 

u(x,0,z) = u^(x,0,z) - l l^g' -g^^^jlKCkj) - EOj)] + 

2tts 


J o 


^ ^ K(kj)l|dx* (26) 

^cJ ■' 


00 

w(x,0,z) = w^(x,0,z) + I ^g'kj [K(kj) - E(kj)] + 


J o 


^ (X-X-) ^ [E(kj) - K Ckj)]]dx' 
^cJ ■’ 


where 


(27) 


kt = 


Sj^ + (x-x')^ + (z - fj')^ 


(28) 


V 2 + k ^ = 1 
^cJ ^ 


36 





and K, E are the complete elliptic integrals o£ the first and second 
kinds, respectively, with modulus k. 

The dynamical condition given by Equation (18) may bo witten in 
the form 


-2 H tU^Cx.O.fj) - 

dx 

1 Vx.O.fj)] 

CjUlc 

1 + 

dxj 

jz 



C29) 


since, in the plane of symmetry 


CYxn) 





1 -.| 

dx 

lVCx,0,£j)l 

1 + 


r 


where 

= I OJ* * U-) 

and 

W = i + w’) 
m 4 

as discussed previously. 

The kinematic condition. Equation (17) , may be ivritten simply as 

df 

(U„ + u(x,0,fj)] ^ - w(x,0,fj) = 0 (30) 

Equations (29) and (30) , with u and w substituted from Equations 
(26) and (27) represent a pair of nonlinear integro-differential equations 


37 



for fj(x) and g(x). The jet location apd slope are knoun at x = 0. The 
function g must display the proper behavior ; that is ^ must vanish at 
infinity as well as having a logarithmic singularity at the origin of 
the jet. Hiis system of equations is to be solved numerically for these 
functions . 


2 . Forms of the Equations for Numerical Solution 
The method adopted here for tl\e solution of the equations dis- 
cussed in the previous section is related to the integral methods of 
17 18 19 

fluid meclianics. * ’ The functions fj and g are represented by 

truncated series with coefficients to be determined in the solution. 

These functions are chosen to give the proper behavior for large x 
and to meet the appropriate boundary conditions. The semi-infinite do- 
main, P ;! x ^ is divided into two regions, a near field 0 ^ x ^ x and 
a far field extending beyond x. The functional representations in the 
near field are of primary concern here. Those in the far field were 
chosen to facilitate integration in that domain; the choice there is not 
particularly critical in the solution since contributions to the inte- 

✓s 

grals from that region are small if x is large enough. Choices for x 
in the computations of 5 to 10 wing chord lengths were found to yield 
satisfactory results. The optimum choice depends on the parameters of 
the problem, such as Cj. This subject will be discussed further in con- 
nection vdth the computational results. 

It may be anticipated that yjy will decrease monotonically with 
increasing x, as it must ultimately approach zero. Moreover, the function 
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dg/dx should display a logarithmic singularity at the origin as suggested 
by the two-dimensional theory of Spence.^ The choice made in this ana- 
lysis is a series of exponentials sucli that 



C 

Ik 


noxp[-S5s) 

I ''J 


I 


+ c. 


InCx/sj) 

M "^x nr/T 


(31) 


SO that 


N 


N-1 


® ' n-\ ‘ ^ 


exp|-If£i - II 




’ x/s , VT r 1 

In — 

• A/Sj + 1 + 1 


InCx/s,) 

— + 2 In 2 

/ X/Sj + 1 

The factor g is a parameter chosen along with x for a particular problem 
to best represent the function over the near field domain. A value of 
g between unity and 10 appeared to be suitable for most of the computa- 
tions attempted. This functional representation is retained in the com- 
putations extending from x = x to infinity as well as in the near field 
region. 

The function fjC^c) was chosen to have the form 


(32) 



(X - 1) 


+ a. 


+ 


(X + 1) 


+ 


% ] 


(S3) 


where x ** vx/Sj and v is a scale factor. Values of v found to yield 
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acceptable results in. the computations attempted were in a range from 
1 to 10. Since the jet passes through the origin, it is required that 


M 


Employing Equation (33) , the slope given by 

V' ^1 ^^2 ^ 

(x + 1)^ (x + 1)^ *” (x + l)^"*'^ 

must agree with that specified in the problem at x = 0. Thus 


(34) 


(35) 


V 




M 

a_j^ - 2 m a^ = - tan (t + a) 
m=l 


(36) 


There remain M free coefficients in the function to be determined in the 
solution. 

Now, while the form of Equation (33) was employed for the jet in 
the near field region x £ x, it was not employed to represent the jet for 
X > X. There the jet sheet was taken simply to have a constant slope, 
matching the value of the near field function at x = x. This represen- 
tation is not precise, of course, but simiDlifies the integrals considerably 
and should not seriously detract from the merit of the solution. Con- 
sidering the inviscid problem, where the sheet would be alloxyed to roll 
up, the resulting trailing vortex pair would have a constant slope, ax^art 
from viscous interactions. 
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Given the above functional forms for g and fj, the final forms of 
the equations discussed in Section 1 may be written as 


df. 


[U„ + u(x,0,fj)] -gi - w(x,0,fj) * p 


(37) 


<(x) = 


2 [U„ + u(x,0,fj) - 

“j 

dx 

w(x,0,fj)] 

Cj c 

r{- 

Ej 

2- 

. 



fN-1 I 

Ln-l " 


n3x 

s. 


In (x/ s j) 

"N . 1^3/2 

(X/Sj ^Ij 


(38) 


uCx.O.z) = u^Cx.O.z) - -J-j. s c j [ ^ k (K - B) 

2ir s, n=l i I 
J o 


dgL* 

^CJ 


(E - K)j dx‘ 


00 

N 

w(x,0,z) = w^Cx,0,z) + 2^ I kj CK - E) + 


(39) 


h 2 1 

.-gr (X - X') (E - K) dx» 

^cJ 


(40) 


and 




t^ere in Equation (29) are now given simply by U,W with the under- 

standing that these are the values computed excluding a suitably small 
region in the vicinity o£ x in the integrals of Equations (39) and (40).* 

B. Method of Solution 

Equations (37) and (38), upon substitution of u, w and k from 
Equations (39) , (40) and (41) comprise a pair of integro-differential 
equations for the functions fj and g. Employing the approximate forms 
for fj and g, these equations become a pair of algebraic equations for 
the M + N coefficients in the solution. With a finite number of terms, 

M + N here, it is not possible, of course, to satisfy tlicse equations at 
all points in the semi-infinite region 0 x < «, Hie equations could 
be satisfied exactly only at ^ (M + N) points at most. Rather than re- 
quiring the equations to be met precisely at the ^ (M + N) points, it is 
advisable to meet the relations in a broader sense over the domain. This 
is accomplished in the analysis by requiring the equations to be satisfied 

A, 

in a least squares sense over the near field domain 0 £ x x, subject 
to values at a larger number of points. 

The mmierical analysis employs a set of control points spanning 

A, 

the near field domain 0 £ x £ x. The number of such points is chosen to 

K fl^ + tl" 

A region small enough such that its contribution to U '= w - ” — is 

negligible. ^ 


42 



be several times the value j 01 + N) such that the least squares method 
will yield smooth functions fj and g, and so that information of the 
necessary resolution within the domain is represented. 

The procedure employed for the solution is an iterative one. It 
may be outlined as follows. 

(a) A trial function k(x) is first diosen. k is employed as 
the primary representation of the jet rather than £j, be- 
cause it is of primary inportance in the dynamical condition 
and it is more sensitive than fj or dfj/dx. Employing k 
at the control points, a first integral obtained numerically 
yields dfj/dx at the control points. These values are 
aiployed to find coefficients in the expression for fj vhich 
yield a least squares fit for dfj/dx, fj and dfj/dx are then 
evaluated from this function at the points employed in ob- 
taining the necessary integrals. 

Cb) These values for f _ and df ,/dx are emiJloyed to calculate the 

J J 

coefficients of the in Equations (39) and (40) at each 
control point. Substituting the results into Equation (38), 
the dynamical condition, the values are obtained so as to 
meet that condition in a least sauares sense over the domain, 
(c) Finally, new '^lues for dfj/dx are obtained at the control 

points, through Equation (37) . The curvature corresponding 
to this new distribution of dfj/dx is then estimated. These 
are compared with the starting values and a new iterate k(x) 
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is generated. When k(x) values obtained in t\ro successive 
iterates agree within a given error bound a solution is 
formally obtained for the problem. 

The scheme described above is not the only one possible, but it ’is stable 
and converges rapidly. It would be possible, for exairple, to employ* the 
linear tangency condition in step (b) rather than the dynamical condition, 
for evaluating the C^. • Then the dynamical condition would yield the new 
K iterate directly. It was found that this method is not always stable. 
Halsey found similar results in a two-dimensional theory. 

IV. NUMERICAL ASPECTS OF THE JET PROBLEM 

The integrals to be evaluated numerically in the analysis require 
careful attention in several respects. They are discussed . together with 
the numerical integration scheme in this chapter. The relations employed 
for evaluating the coefficients in step Cb) of the iteration scheme 
are then summarized. A description of a method of evaluation for the 
starting values and the new jet function is then given folloxved by a 
discussion of the procedure adopted for generating a new iterate k(x) . 

A. Evaluation of the Integrals 

The complexity of the integi*als involved in the equations for 
the velocity field requires that they be computed numerically. Three 
difficulties arise in this connection; 

1. singularity of teims in the integrand at x' = x, the con- 
trol point 
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2. infinite v^per bound 

3. singularity at the origin caused by the logaritluni.c term in 
the. vortex strength expansion 

1. Singularity at the Control Point 
At the point x' = x, the modulus of the comi)lete elliptic integrals 
(ICj) becomes unity and the elliptic integral of the first kind becomes 
infinite. This presents the problem of niunerically integrating through 
a point in vdiich some of the terms of the integrand arc singular. Of 
course, the integral is finite. It is possible to circumvent this 
difficulty by applying the identities 



Inserting Equations (42) through (44) into the integrals under considera- 
tion results in the expressions 
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uCx.O.fj) - u^(x,0,fj) hrj^^^nj 


«n' «j’ >'cJ O' - B) . 


SjTF 


1 + 


X - X 


'II 


■»V V 1 - 

X - X' ,2 - t /.J 

imr -‘ 


dx' 


X - X’ 


/■4 r% 


and 

wCx,0,fj) » w^(x,0,fj) + 




8n’ ’^cJ O' - E) 

■^•’Fisir 


■iSn' 1 ® W 

T3?-ig 


rr^ 


1 + 


J 2 


X - x' 




X - x' 


(46) 


Since 


lim ^ 
k-^-l 


(K - E) -»■ .0 



the e:!qpressions vrithin the braces of Equations (45) and (46) are con- 
tinuous and non-singular in a small neighborhood of the control point 
(x' = X ). These teims- are plotted in Figure 4 as f and f 2 » respectively. 
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The associated integrals may be written as Cauchy principal values. 


f 


X“C 


1 


f.Cx*) 

-i r 

V - x* 


; i = 1,2 


in a neighborhood of x* » x. It is then possible to expand f-^(x') in 
powers of (x - x’) within this neighborhood and to perform the integration 
analytically. By choosing and to be small enough, this expansion 
may be truncated to fourth order in (x - x') without loss of accuracy (refer 
to Appendix F). In any event the contribution of this element to the 
total integral is extremely small. 


2. Infinite Upper Limit 

The problem of numerically integrating to infinity may be removed 
by employing a change in variables. A natural candidate for this task 
is the modulus of the ellpitic integrals kj. From Equation (28) it can 
be seen that 

s,^ dkr 

dx' = — cff-T 

(X ^ X') + (z - fj*) 

A 

By en^jloying Equation (47) an integral over the range x x' £ «» may be 
transformed into an integral with respect to kj over the internal kj £ 
kj ^ 0 vdiere kj = ICj(x) . For exanple 
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/ (K-E) * \ (E-k^^ X)] dx- 

* l=cJ J 


V 

■4 


’ 


K-E , 

"-V S-kc/ Kl 


oc 

x-x' 

dx' 

^2 dx' 

X-X‘ . 2 , 2 

*^cJ J 

1 + , 

rxi 



W 

dx' 

^ x-x*^| 



By noticing that 


lim K - E _ iM ^ ^ ^ ir 

1(^0 k>0 jj.2 i T 


and 

x-»« X - x' ~ x*o» dx' 

it can be seen that the resulting integrand is zero at kj = 0, The slope 
dfj/dx is taken to be a constant here, in accord id.th the approximate 
representation of the jet in this region as discussed in section A2 of 
the previous chapter. The resulting integrand for the x (z) velocity 
component is plotted as f^^ 0^3) ^ Figure 5. 

3. Singularity at the Origin 

The final difficulty encountered within the integrands of interest 
involves the logarithmic singularity at the origin due to the term dg^^/dx. 
Consider only the integral involving this term in a small S neighborhood 
about the origin. 


4a 



C«) 


f - fj 1 B - ''rf K tel 

J T5r X - X’ ^ — Tz-Tjp- x-x 


rhe group of terms 


x-x' 


F" 



7 



(plotted as in Figure 6) may be given their value at a mean point 
and removed from tlie integral to provide the resulting appi'oxunation 
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A similar relation was developed for the z velocity component (the term to 
be evaluated at the intermediate value is plotted as €2 Figure 6) . 

TTie integrands of Equations (39) and (40) have been plotted as 
and £2 respectively in Figure 7 for the remaining intervals. 
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4. Numerical Integration Method 
In accord vdth the factors mentioned above, the integrals were 
subdivided into at most the sum of five integrals over parts of the do 
main of interest. 


p < X* <6 
6 £ X' ^ X - Cj 
X - < X* £ X + £2 

X + £2 £ X* £ X 

yv 

X < x' < «» 


where 6, £^^, and £2 are small values chosen on the basis of integration 
accuracy. Those integrands not involving the logarithmic singularity 
at X = 0 were integrated directly from 0 to x - 

Now, in order to represent the integrals well for a range of jet 
geometries, it was decided to choose the main control points at nineteen 
equally spaced locations along the jet (A5 = constant). An additional 
control point was placed at A?/2, within the initial region of rapid 
variation of the jet geometry, resulting in a total of twenty control 
points in the domain 0 £ x' £x. Step sizes for the integrations were 
taken to be ten times smaller, resulting in a total of 201 points in 
the numerical grid. These values were found to give satisfactory accuracy. 
Numerical integration of the functions with nonsingular integrands was 
perfonned by means of a compound Sin^son's rule.^^ 
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J fCx) dx - y [fj + 4 Cfi + fj + ... ♦ fjn.p ♦ 2 Cfj * «4 ♦ ... + 

V2) 

v^re 

fi = f (xp 


h * (Xj^^.2 ~ x^) = constant 
and the remainder is given by 


\=- 


nh^ 

T!r 


d^fK) 

dx^ 


< 5 < X. 


2n 


Hie indices refer to the grid points within the integration interval 
of concern. Numerical algorithms for the complete elliptic integrals 
are described in /^pendix B. 

B. Evaluation of 


Upon evaluating the integrals numerically at each of the control 
points, one has 


N 


u (Xi) = (Xj) ♦ C„ I„ (Xj) 


and 


N 


w (xp = (X. ) + C„ J„ CXi) 
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vdiere denotes the control point and (x^) are now known 

coefficients. These expressions may be substituted into the dynamical 
condition, Equation (38) to yield 
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(50) 


This system of twenty equations in N unknowns (the N values) is an 
overdetermined system (N was chosen to be 6 in the computations). Tlie 
values are evaluated such that Equations (50) are satisfied in a least 
squares sense over the domain. 

The least squares solution may be obtained by minimizing the sum 
of the squares of the residuals 
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(53) 


and 


52 




Substituting the expression for Q into this relation yields 



which may be solved for the C values. 

' n 

Equations (56) comprise a nonlinear system of algebraic equations 

for the N values of the coefficients A solution may be obtained by 

21 22 

employing a variation of the secant metliod ’ . Ihis is an iterative 

procedure starting from a trial set of values ^ for most such 

nonlinear problems, the trial values for stai'ting tlio solution must be 
reasonably close to tlie solution for rapid convergence of the method. 
Once the computations are undeiway, preceding iterates serve this pur- 
pose well enough. In a new computation, liowever, it may be necessaiy 
to provide initial values v/ithout tlie guidance of prior solutions. One 
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means for accomplishing this is simply to employ tlie linear tangency 

condition, Equation (37), solving for the .C^- values enn^loying the 

current fj function. These values of may then be used to calculate 

the velocity, components , u and w. Since tlie velocity components change 

n 

less rapidly than the vortex strength , a linearized dynamic condition 

may now be employed (with only y as an unknoun) to obtain an improved 

set of values {C }. This last set of C values can then be employed as 
n n 

a trial set in the solution of the nonlinear dynamical condition. 
Experience with this procedure has been satisfactory. 

C. Evaluation of the Function fj 

The function fj is determined iteratively in the solution pro- 
cedure. An initial trial value must be chosen for a new problem. Sub- 
sequent values are generated by the program in the iteration scheme. 

The procedure, outlined in the previous discussion, amounts to (a) 
starting with a trial function fj, (b) evaluating g(x^) given fj so 
that the dynanical condition is met, and (c) employing the tangency con' 
dition to guide the choice of a new iterate for the function fj. Both 
fj and g must be such as to satisfy both the tangency and dynamical con' 
ditions in the solution. 

Upon initiating a new comixjtation, it has been found convenient 
to start with a function f_ derived from the t^vo- dimensional solution 

k) 

of Spence.^ The appropriate equations for this purpose are listed in 
^pendix C, Wliile this method is best for small Cj, a, and t values. 
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it is generally preferable to employ solutions obtained in the present 
analysis for similar parametric ranges, v;hen they are available, if Cj, 
o, T are large. More rapid convergence can then be adiieved, as would 
be expected. 

Tlie function fj giving the position of the jet in the analysis 
is not itself sensitive enough to be employed either for the estimation 
of errors involved or for iterative purposes. Since both the first 
derivative and the curvature involving the second derivative arise in 
the theory they must be represented with sufficient accuracy and regular- 
ity in the numerical scheme. Tlierofore, the curvature, being most 
sensitive, was employed both as a basis for the choice of new iterates, 
and as the factor on \diich accuracy was judged. 

The computations start with a set of values for at the 

twenty control points. This function is numerically integrated to yield 
values for dfj/dx at the control points. Hie method employed for inte- 
gration of the curvature may be found in Appendix D. The coefficients in 


the expansion of Equation (34) for f j are determined so as to provide a 
least squares fit to the slope, dfj/dx, at the control poiiits. Thus one 
has25>24 


J M M 

2 a. [ S G.(x.) G^Cx.)] = Z F(xJ G,.(x.) ; k= 1,...,J (57) 

j=l J i-1 J ^ ^ ^ i=l 1 K X 


where a^ (j « 1, J) are the coefficients in the expansion for f 
(cf. Equation 33). and 
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G^(x-)=j 1- 


(58) 


(Xj + ] 


F(x.) = - 


dfj df j (x = 0) 


dx 


( 59 ) 


Upon deteimiiiing tlie coefficients a. in the representation of 

3 d 

it is necessary to determine positions of the new grid ppints Xj^ such 
that they are equally spaced along the jet. Equal spacing of the grid 
points is desirable for better representation of the fiuictions involved 
in the numerical work as mentioned in Section A4, Tliis amounts to 
numerical integration at the expression for d?/dx 

A A 

o o 


dividing the interval K iJ^to suitable equal intervals (200 here), and 
establishing a correspondence between and Xj^ for each grid point. The 
relations involved and the integration scheme anployed are described in 
Appendix E. 

Ihe computations employing fj proceed as discussed in the previous 
section to an evaliiation of {C^} satisfying the dynamical condition. The 
tangency condition remains to be met. This relation 

♦ u(x,0,fj)] -JJJ - w(x,0,fj) - 0 

may now be employed directly to evaluate dfj/dx at each control point. 
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D. Generation of a New Iterate 


As mentioned in Section C above, it is advantageous to employ k(x) 
as the basis of the iteration scheme so that tlie key aspects of the jet 
with respect to the mechanics of the problem may be represented properly. 
Values of k employed at first in the computations arc smoothly varjdng 
in the domain. New iterates are also to be required to vary smootlily. 

The problem here, of course, is tliat of generating smooth clianges 
for K based on information at the discrete control points regarding dfj/dx. 
From the definition of tc, a small clumge may be related to changes in 
dfj/dx and d^fj/dx^ as follows 



A procedure for obtaining a set of 6k values at the control points which 
proved effective in the iteration procedure was simply to ignore the 
first term (it could be anticipated to be small), and to employ the re- 
lation 



wlvere k^^ represents the value of k employed in starting the preceding 
iteration and 6(dfydx). represents the difference between the computed 
values of dfj/dx tmd those values employed in starting the precaling 
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interation. Hie iteraticn procedure converged when the rtnxv iterates were 
chosen according to tlie relation 

where e is a damping factor cliosen to be small initially and to increase 
to unity as the solution was approached. Initial values for e in a new 
problem might be on the order of 0.2, 

A measure of the error in the result may be obtained by comparing 
these new values of at each control point with those employed upon 
starting the previous iteration, A sum of the absolute value of the re- 
lative errors between these two values was employed in the computations for 
this purpose. Hie existing method employed of value of 0.6 for this sum 
(or the average relative error at each point was required to be below 0,03), 
If this criterion was not met, the neiv values for the curvature were 
employed to initiate another iteration, 

V. ANALYSIS OF THE FLOWFIELD OF THE AIRFOIL 

There are two aspects of the aerodynamic field of the airfoil 
sntering the theory. The first is the need for some initial estimate 
Df the influence of the' flowfield induced by the wing on the jet. This 
nay be obtained by means of one of the sinpler linear wing theories 
for exanple. The second is the solution of the flowfield about the 
airfoil subject to the full influence of the disturbances introduced 
by the jet. The first is treated here assuming the airfoil to be 
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represented as a lifting line vdth the appropriate trailing vortex dis- 
tribution. An elliptic loading is assmed for tliis purpose with a total 
lift estimated to represent the final solution (at least somewhat larger 
tiian tliat for the wing alone) . The second is accomplished by numerical 
methods in an adaptation of lifting surface theory. 


A. Initial Estimate of the Flowfield of the 
Airfoil for Use in the Jet Solution 

For an elliptically loaded wing operating with a lift coefficient 
Cj^, the distribution of the vortex strength along the lifting liiie may 
be expressed as^"^ 

■• = r,yrT(^ C61) 


where 


A z L o® 0 


The vorticity associated witn the bound vortex may be placed at the 
quarter cliord point (x = Xj^) of the airfoil^^ and the trailing vorticity 
may be regarded as lying in the z = 0 plane (refer to Figure 8), 

The velocity components u^ and w^ induced by tliis vorticity dis- 
tribution may be determined from the relations 

^A 


u. = 
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(X - 
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(62) 
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dy' (63) 
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Values of and at y = 0 are required in the jet solution. These may 
be obtained by carrying out the integrations of Equations (62) and (63) 
subject to the vortex strength distribution of Equation (61) to yield 
(cf . Appendix A) 


u^(x,Q,z) 
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z k 
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W^(X,0,Z) . - 
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(65) 


vhere K and E are the first and second complete elliptic integrals of 
modulus k^ and II is the tliird complete elliptic integral, whidi is also 
a function of the parameter And 





(x - xj^ + s.^ + 
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B. Solution of the Flowfield about the Airfoil Subject 
to the Interference Field of the Jet 

Given the solution of the jet problem as described in Chapters III 
and IV, subject to the estimated interaction of the airfoil of Section 
A above, it is next required to obtain a more precise solution for the 
flowfield about the wing. This is formally given by the relations 
(recall the discussion of Section II B-3) 

Jir II <Ka* - V, (J) (66) 

®A 

subject to the kinematic (tangency) condition 

af. af. 

at points on the airfoil whose surface is given by tlie equation 

- 2=0 

It should be pointed out that unlike linear lifting surface theories 
the streamwise disturbance velocity is no longer negligible compared to 
the freestream velocity. 

The velocity components Uj and Wj at points on the surface of the 
airfoil may be evaluated by integration over the vortex strength distribu- 
tion of the jet 
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■ w II 




( 68 ) 


Only airfoils with little or no dihedral (i.e. 3f^/3y is negligible) are 
considered here. In addition it is assumed that the spanwise velocity 
component is small conpared to the do;vnwash and the streaim\dse velocity 
conponents. Therefore, the equation necessary for the evaluation of this 
COTponent is dropped from the ensuing treatment. In Eqiiation (68), the 
integrations xvith respect to y' may be obtained analytically (refer to 
Appendix A) . The remaining integrations must be carried out numerically 
in a manner already described in Chapter III, 

It remains to arrange the surface integral of Equation (66) , which 
gives the contribution of the vorticity associated with the airfoil to 
the velocity field, in a form suitable for solution of the airfoil problem. 
Ihe surface over which this integral is evaluated includes the surface 
used to represent the airfoil (the projection of the wing onto the z - 0 
plane) as well as the surface containing the vorticity which is waslied 
do;>nstream, the wake vorticity. The location of the wake vorticity is 
not known, a priori . For a relatively lightly loaded airfoil the contri- 
bution of the vorticity created by the airfoil may be considered negligible 
ccmpared to the vorticity of the jet in this wake region. However, in 
order to simplify the resulting equations the vorticity shed by the airfoil 
will be retained. It is assumed that the location of this vorticity may 
be taken to lie in the z * 0 plane. This assumption is made in most air- 
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foil theories. It is therefore possible to utilize some of the mathe- 
matical aspects of those theories* 

Bnploying the definition of a vortex sheet the vortex strength 
of the airfoil may be presented by the equation 

= -Av e^ + Au Oy. (69) 


vdiere 


Au = u - u 


Defining a doublet strength as 


= « - $ 


♦ being the velocity potential. Equation (69) may be vnritten as 


Ya = - 


3y 




3x ®y 


(70) 


Equation (70) along with the assumptions discussed in the preceeding para- 
graph may now be employed to provide the resulting equations for the 
induced velocity field created by tlie airfoil 
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Solution of the airfoil problem requires these equations be applied to the 
z = 0 plane. It may be seen from Equation (71) that = 0 for this plane; 
therefore, only Eqpaation (72) remains to be evaluated. 

Equation (72) as it stands contains a singularity with respect to 
both integrals at the control point (x,y) ; therefore, a more suitable 

25 

kernel will be sought. Following the procedure of Robinson and Laurmann 

"d' -fz (eI 

3z 


this equation may be witten as 

CO 

" W / \ 


A.e. -s. 


Further manipulations of Equation (73) provide the final working equation 

26 

for the downwash induced by the airfoil in the fonn ^ 
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(74) 


for an arbitrary field point (x, y, z) . In the limit z 0 the second 
term of Equation (74) vanishes and the first is given by the Cauchy prin- 
cipal value 
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Wj^(x,y,0) ■ - Jj. 



3x' 


(1 + 



- y 
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By defining 


f(y’) 
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(75) 


employing integration by parts, and expanding £(y ± e) in a Taylor Series, 


Equation (75) may be obtained in the final form' 
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w^(x.y,0) = I f ^ Cl * dx. 

_s ly - y J ji.e. 


(76) 


where the symbol ^ represents tlie Mangier principal value. The Mangier 
principle value is defined by the equation 


“S. f J -s ^ ^ y+e 
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Equation (76), found in many of the existing lifting surface theories, is 
singular with respect to only the outer integral. The integral with re- 
spect to x' possess only a discontinuity, since for y =* y', 

lim r, X - X*, J2 : X > X* 
x-»-x' ^ R |0 : X < X* 
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The integral equation. Equation (66) , \d.th given by Equation 

(76), together with the kinematic condition. Equation (67), now are in a 

form suitable for matliematical solution. Given a solution (x,y) on 

the surface of the airfoil, the floivfield may be computed at a point of 

interest. There are several methods which may be employed to obtain this 

desired doublet distribution. The method adopted here is related to the 

27 

finite element theories 

The pressiire difference across the airfoil may be obtained from 
Bernoulli's Equation. 

ip - I p C-2(JI„ + Uj) - 2 Vj iv^ t . 

2Wj 

As a consequence of the fact that the singularities used to represent the 
airfoil have been placed in a plane, the resiilting velocity potential for 
a cambered wing at incidence to the freestream is antisymmetric with 
respect to that plane (z = 0 plane). Thus, 



Therefore, the pressure difference across the airfoil may be expressed m 
the form 


Ap = -P + Uj) Au^ 


(77) 
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vfliere the assumption that the spanwise velocity component is negligible 
con^ared to the streamiidse component has been employed in this last 
equation. 

Utilizing ■^e definition of the doublet strength along with Equation 
(77), a relation for the change in the pressure coefficient across the. 
airfoil may be determined 


ACp = - ^ ^ 


U' 


ax 


(78) 


00 


Substituting Equation (78) into Equation (76) provides the result- 


Consider the airfoil to be canposed of a number of finite, recLaugu- 

lar elements. In a manner similar to finite element theories, it is 

assumed that the pressure difference across each element is constant. In 

addition, it will be assumed further that the velocity field induced by 

tlie jet is also a constant on each element. The downwash due to the 
th 

element at the i control point can Lheu be written (refer to Figure 9) 
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( 79 ) 


The resulting surface integral in Equation (79) is only a function of the 
geometry of the element. The evaluation of this integral is treated in 
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Appendix A where proper care has been taken for the situation in which k 
i. Defining 





the total downivath at the i™ element » for the case in which the wing 
planfoim is divid^ into 2N elemoits, may be written 
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Because of geometric symmetry the pressure difference and velocity field 
induced by the jet must be symmetric about the x~axis, therefore 
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(81) 


th. 

where the k element lies on the right semi-span wing section (y > 0) and 
has as its reflection about the x-axis the k + element (refer to 
Figure 9) . It may be shown further that 


'k+N,i ' - ='ak' ^bk’ ■ ^ak- - ^bk’ 


The unknown pressure difference for eucli element (AC^j^) may now 
be obtained by applying the kinematic condition at the N control points, 
which have been taken to be located at 
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Xi = X(,i + 0.1 (x^. - Xj,p 


O.S 

Inserting Equation (81) into the kinematic condition results in a set 
of N simultaneous, linear equations for the N unkno^vn pressure differences 
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C. Calculation of the Forces on the Airfoil 


Once the pressure distribution on the airfoil has been obtained the 
pressure forces may be computed directly. 


where 




Ap' cos ijf' dS^' 


Ap' sin \p' dS^' 





(83) 


(84) 


(For the situation imder consideration e a). In the above equation 
Op) refers to the lift (drag) force experienced by the airfoil due to 
the pressure difference across it (as opposed to tlie force created by 


69 



tfie reaction to the jet). It should be mentioned that Equations (83) and 
(84) do not include the leading edge suction force which provides only a 
negligible contribution to the lift. Since the flow is considered inviscid, 
the above equations also exclude skin friction effects. Written in terms 
of non-dimensional variables, Equation (83) has tlie form 




cos T|»' AC * 
1 


where is the surface area of the airfoil. In terms of finite elements 
covering the surface. 



1 2N 
— T. 

k=l 


cos 


C85) 


where 

\ ' ^^bk " ^ak^ ’ ~ ^ak^ 

is the area of the k element. It should be noticed that the assun 5 >tion- 
that the angle ip is constant for an element has been employed in Equation 
(85) . Such practice is inherent to the definition of a finite element in 

on 

theories such as Woodward 's.“' Utilizing the symmetry conditions. Equation 
(85) may finally be written in the form 

^LP “ " ^ 1 ^k 
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Similarly, the drag coefficient resulting from the pressure difference 
across the wing may be obtained from the equation 




2 ^ 
k=l 


sin 


k \ 


( 87 ) 


The total lift and drag coefficients for the airfoil may now be 
obtained by including the reaction to the jet in Equations (86) and (87) 








J Cj c sin (t + a) dy 
■®J 

I Cj c cos (x + a) dy 


( 88 ) 
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D . Calculation of the Velocity Field in the 
Vicinity of a Tail Surface 

The dovmwash at a point in the vicinity of a tail or control sur-- 
face downstream of the wing may be immediately obtained from the equations 
of the previous sections. For simplicity, a point located on the x-axis 
is anployed for this purpose. Tlie velocity at such a point is written 



Again the first surface integral will be referred to as u^(w^) and the 
second as Uj(Wj), 
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1 • The Velocity Field Tiicl uced by the Jet 


The necessary equations for tJie cletennination of the velocity field 
induced by the jet sheet at the point under consideration may be obtained 
from tlie formulation of Qiapter III. Because of symmetry the sprinwise 
velocity component is zero. Furthermore, 
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Proper account is taken for the infinite upper limit and the singularity 
at the origin in exactly the same manner discussed in Chapter IV. 

2. The Velocity Field Induced by the Airfoil 

The velocity induced by the airfoil at the point (x,0,0) may be 
obtained by employing the results of Section B. 
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U^(x,0,0) = 0 


w^Cx,0,0) 


v^(x,0,0) =* 0 



(90) 


Bnploying finite elements Equation (90) may be witten in the form 


w^Cx.0,0) = - 


N 


AC, 




“ “jk-' 


where Ij^ may be found in Appendix A. 

VI. COMPUTATIONAL RESULTS AND COMPARISON WITH EXPERIMENTAL 

DATA AND OTHER THEORIES 


A. Choice of Parameters for the Computations 

The airfoil and jet geometries as well as values for the angle of 
attack of the wing, the jet deflection angle, and the jet momentum coef- 
ficient were chosen largely in accord with data available in the literature 
for the purpose of comparison. Computed results were obtained for a 
rectangular airfoil ^\dth an aspect ratio (.AR) of 6.8 equipped with a full 
span jet flap. Choices for the remaining parameters along with the com- 
puted results for the lift coefficient, the induced drag coefficient, and 
the downwash angle (0) and pressure coefficient (Cp) three chord lengths 
downstream of the airfoil are listed in Table 1, 

A first solution was obtained for small values of a, t, and Cj by 
employing initial values for the jet geometry obtained from the tivo- 
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dimensional theory of Spence^ (refer to ^pendix C) . Subsequent soluticms 
were obtained using initial values obtained from the results of an existing 
solution which had input parameters (a, t, and Cj) close to those under 
consideration. 

Choice of the number of grid points employed to determine the neces- 
sary integrals and the number of control points on the airfoil and the jet 
was made on the basis of providing the highest accuracy for reasonable ccmt- 
putational times. By using 201 grid points, 20 control points on the jet, 
and 32 elements on the wing, execution times of approximately 26 seconds 
per iteration of the jet problon and 23 seconds to obtain the resulting 
forces on the wing and induced downwash at a downstream location were 
obtained. All confutations were made on a CDC 6500. Hie factors $ and v 
(refer to Equations (32) and (33)) were each given a value of 10, for the 
computations. The resulting values for and a^^ fell in the ranges 

0.05 < IC I < 6 
* n* 

10 < Ja^^l < 160 

The nearfield calculations were obtained in the region extending from the 
trailing edge of tlie aiffoil to a point (x) approximately 5 chord lengths 
downstream. It will later be observed that this region contained the 
primary contribution to the vortex strength, as it should. 

A solution to the jet problem was assumed to be determined once the 
sum of the magnitudes of the relative errors between the curvature values 
obtained from two consecutive iterations was below 0.6 (resulting in an 
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average relative error of 3 percent) . However, for larger values of Cj 
and a it was found that the relative errors of the curvature nt each con- 
trol point were approximately constant prior to satisdfy^ing the above con- 
vergence criterion. Consequently, when tliese values \trere scaled to pro- 
vide proper integrals (refer to Appendix D) the resulting changes from tlie 
previous iteration was small. In such cases convergence became slow when 
the average relative error became less than about 7 percent. An investi- 
gation into this dilemma showed that the application of a decreasing fimc- 
tional form for the danping factor (e) was fruitless. However, by assuming 
that the value of the jet slope at infinity (refer to Appendix D) was half 

/N. 

its value at x (a value of zero was employed in the original conputations 
for convenience) , it was possible to reduce the average value of the 
relative errors to 2 percent. Since the resulting value for the relative 
error of the lift .coefficient was only 3.5 percent, however, the above 
error criterion was reduced to a value of 1.3 for these cases. Iterations 
of the entire problem were terminated once the relative error between two 
consecutive values of the lift coefficient was below 0.025. Ihe. above 
convergence criteria required approximately 3 iterations of the entire 
scheme and an average of 6 iterations of the jet scheme. For the initial 
values, however, as many as 10 iterations were required to converge on a 
jet solution. 
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B. Results 


The variations o£ the lift coefficient with changes in the jet 
momentum coefficient are provided in Figure 10. Comparison with the linear 
theories of Nhskell and Tokuda^ are favorable. For small values of Cj 
(linear region) the disagreement between the theories was less than 10 
percent. However, the values of computed by the present theory for 
large values of C (nonlinear region) are evidently higher than wliat would 
be predicted by extending the cun^es obtained fi'om the linear theories. 
Experimental data^^ are also available for the linear region. It can bo 
seen in Figure 10 tliat agreement between tlie present tlieory and the exi^eri- 
mental data is also quite acceptable in that regime. It should be mentioned 
that the results of Maskell have been increased by an appropriate factor 
to account for the thickness of the airfoil employed in Alexander's^^ experi- 
ments. The factor involved is discussed in greater detail by Lissaman. 

It can be seen in Figure 11 that the resulting variations of Cj^ 
with changes in a provided by the present ( nonlinear ) theory lies slightly 
below the experimental data and the curves obtained from tlie linear 
theories. However, it has been observed that linear theories in general 
tend to overestimate the lift, the difference increasing with larger values 
of a and t. The linear tlieories employed a lifting-line approximation 
for the airfoil. It might be expected that the utilization of a lifting- 
surface technique vrould lower these values in a manner sbiiilar t'o that for 
conventional wings. It has been suggested that the linear theories pro- 
vide favorable agreement with the experimental data because of the "fortui- 
tous cancellations" of higher order effects.^ The results of Lissaman’s 
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theory are also depicted in Figure 11. They tend to lie even lower than 
those o£ the present theory. 

The results for the induced drag coefficient are presented in 
Figures 12 and 13. As already discussed in the previous chapter, the in- 
duced drag Cor thrust for the conditions encountered) resulting from 
the present theory only accounts for the reaction to the jet and the 
pressure distribution on the airfoil. Contributions to the drag attributed 
to form drag and to viscous effects, such as skin-friction or jet entrain- 
ment, are not represented in this analysis. Drag coefficient values for 
a “ 0® as well as values obtained from Figure 12 for a 0® are presented 
and compared to Williams data in Figure 13. The agreement is again 
quite acceptable. 

The distribution of the pressure difference across the wing is 

shown in Figures 14 through 19 for three different cases. Since only four 

elements were employed in each of the spanwise and chordwise directions 

only basic trends can be obtained from these figures. For more refined 

results it would be necessary to employ more elements at the cost of 

computing time. As expected, a larger than normal pressure difference is 

encountered in the vicinity of the trailing edge of the airfoil. For this 

reason a wing equipped with a jet flap experiences large nose-down pitching 

moments. The spanwise distribution, however, is not very different from 

10 28 

conventional wings. This has also been observed in experiments. * 

The locations of the jet sheets in the plane of symmetry are shorn 
in Figures 20 and 21 for different values of Cj and a. Obviously, the 
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jet lies well below the z = 0 plane. For this reason downwash results ob- 
tained from linear theories are open to serious question. 

The vortex distributions provided in Figures 22 and 23 display the 
proper logarithmic singularity at the origin. It is also apparent from 
these figures, that the control points were located in a region which 
contains most of the vorticity. The distributions indicate that the vortex 
strength becomes slightly negative beyond x = Sj. It is felt that this is 
a result of errors in the truncated series approximation enq)loyed for the 
vortex strength; it is not likely to be real. 

As a result of the large nose-down pitching moment experienced by 
an aircraft equipped with a jet flap, an exceptionally large amount of 
trim is required. Consequently, a control surface plays an important role 
in such configurations. It is important, therefore, to know tlie flowficld 
with some degree of accuracy in regions where the tail surfaces would be 
located. Downwash results obtained from the present tlieory are shorn in 
Figures 24 through 27 in tlie form of tlie downwash angle (0 = tan ^ (~U/^V)) 
and momentum or pressure information (Cp) . The results were obtained at a 
representative point along the x-axis 3 chord lengths downstream of the 
trailing edge of tire airfoil (x = 0.88 Sj). 

VII. CONCLUSIONS 

The results of the present theory show ample agreanent with the 
available experimental data and compare well rvith existing linear theories 
in the operating regime where they apply. For the lift coefficient, both 
the variation with Cj and a found in the theory were satisfactory. The 
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values o£ Cj^ con5>uted were slightly below the experimental data by an 
increment which could easily be accounted for by considering the airfoil 
thickness effect both in magnitude and trend. The induced drag (thrust) 
coefficient also agreed for the conditions considered. Regrettably, there 
are no experimental data available for the position of the jet sheet for 
a v.ung equipped v/ith a jet flap for larger values of Cj, a, or t. This 
situation should be remedied in experimental research. In addition, tliis 
theor)'^ extends values for the force system of the airfoil over a wide 
range of values for Cj and a, beyond that for which the linear theories apply. 
Moreover, the downwash computations are likely to be a substantial iiiqjrove- 
ment over results obtained from the linear theories. 

The existing theory easily may be extended to include the increase 
in the lift due to the thickness of the airfoil and the jet. TTiis may 
be accomplished by employing a proper source distribution over the vortex 
sheets used to represent tlie actual wing and jet. Furthermore, it has 
been observed that entrainment effects are capable of providing substan- 
tial gains in the lift coefficient under some circumstances. This is 
especially true for the situation in which the jet is exhausted over a 
small flap. Inclusion of entrainment phenomena into the present theory 
is also possible. This may be acconplished by locating a proper sink 
distribution over the surface used to represent the jet. 

It should be mentioned that this theory is not limited to rectan- 
gular airfoils equipped with full span jet flaps for vdiich data was ob- 
tained in the previous chapter. However, as it exists no account can be 
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made for spanv/ise variations in Cj and x or a wing with dihedral. More 

general jet properties may be included by employing a suitable functional 

form for the spanwise variations of the vortex strength (gj .* Furthemnore, 

the present theory was developed in a manner that any wing theory might be 

employed. Therefore, more general airfoil geometries may be handled by 

27 

utilizing more sophisticated wing theories such as Woodward's. Finally, 
it should be noted that the ground effects problem may be treated by 
employing an imaging technique. 

The present theory contains many parameters such as the number of 
control points, the grid point spacing, the factors 3 and v, etc. whose 
values must be specified ab initio . An investigation of the effects 
changes in each of these parameters has on tlie solution procedure and re- 
sults should be undertaken. Ultimately, it would be desirable to develop 
an algorithm v;hich woiOd automatically adjust the values of these constants 
depending on the values of x, a, and Cj. Furthermore, a study of different 
functional forms for the damping factor (e) which would achieve convergence 
at a faster rate should also be considered. 


but not without some additional mathematical complexity in the integrals. 
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APPENDIX A 


Evaluation of the Integrals 

Hie theory developed within the main text contains several integrals 
o£ standard form. A brief outline of the integration procedures and a 
summary of the results are provided herein. 


A. Integrals Inherent to the Horseshoe Vortex System 
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1 ) 




dn 
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Referring' to Equation (213.06) of Byrd and Friedman for the case m “ 1 
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B(i 5 >loying the equations found on pages 45 and 49 of Petit Bois 

I, - [sin*^ (n) £ ^ -^l ^ ? [l ] 


(AS) 


3) 


X 

-=l 


nidn 


0 (n^ + q) 7(1 - n^) (n^ + c^) 


(A6) 


Consider Equation (213.12) o£ Byrd and Friedman 
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where k satisfies Equation (A3). 


B. integrals Inherent to the Jet Problem 
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Taking advantage of the symmetry of the integrand and employing the trans- 
formation 
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such that 
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then 
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Both of these integrals are elliptic and may be found in Byrd and Friedman 
(Equations (235.00) and (235.01)). 


- 2k [K(k) - E(k)) (A9) 

idiere k satisfies Equation (A3) 
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Consider the transformation 
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such. that. 




The original integral may then be writtoi 
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From Equations (218.10) and (318.02) of Byrd and Friedman' 
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where k satisfies Equation (A3) 
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Bnploying Equation (259.05) of Byrd and Friedman' 
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where 
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Furthermore, it can be shown that 
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\diere r^ satisfies Equation CM3) above. Bi^loying the identity 
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Equation (A19) may be witten 


. .--i 

n 3 y J 


n dn 


,{ yci - nui + nJln - c)(n - c-^) 


Referring to Equation (259.03) of Byrd and Friedman 
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when A, B, g, and k satisfy Equations (A15) through (A18) and 


(A20) 
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A - B 
A + B 


« - A + B 
“2 ■ A - B 


(A21) 

(A22) 


By utilizing the "Addition Equations" (Equation (10) of Byrd and 
30 

Friedman) it can be shown that 

n (», k) - 2 n (Y^, k) 


Employing the additional relations 
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the resulting foim for Equation (A20) may be obtained 
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It can be shorn, with some difficulty that 
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C. Integrals Inherent to the Finite Element Teclmique 


Ic H i , f (1 + 5^-E) a (A28) 

\diere 9 

r^ = (X - + Cy - n)^ 

Defining the inner integral as J 

J = C + f - - - ■■ ^ dC 

7(x - + (y - n)^ 

Using the transformation p = (x - C) » the resulting integral may be found 
in Dwight^^ 

J = 5 - VCx - if + (y - n)^ (A29) 
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in vdiich 5 is given the value of the limits. 

The remaining integral must be considered over two intervals, one 
craitaining the point y = n and one not containing this point. The latter 
situation provides no difficulty 
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The first integral can be found immediately 
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The second integral may be determined by employing the transformation 
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The final result may be written 


h “ ^ ' ^ - .*,0^ ' * 1” [Cy - n) ♦ V(x - c)^ ♦ (y - n)^ ] (ASO) 

If the interval of integration includes the point y * n one must 
handle the resulting singular integrand. This may be accomplished by 
employing the Mangier Principal Value which is defined as 

4 fCn) dn ^ lim f f f(n) dn . f f(ri) dn . 2 f(y) 1 

I Cy - I i (y - y ieCy - J 

It may be shown that the resulting limit is identical to the result given 
in Equation CA30). 

By substituting arbitrary constant limits into Equation (A28) , the 
functional form en^loyed in Chapter V may be obtained 

d b 
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+ In y - c + (X - a)^ 

y - c + - c>^ + (x - b)^ 

- 3^ y-d. Jcy-df* (x ^af 

y - d + 7(y - d)^ + (X - b)^ 

Although the integral anployed for the doimwash calculations is a 
special case of Equation (A28) use was made of the fact that the integrand 
was symmetric. The integrations were therefore performed over an element 
and its reflection about the x-axis. 
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vdiere use has been made of the identity 
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For elements lying along tlie x-axis (c = 0) the Mangier Principal 
Value"” must be obtained. The resulting form is 

d b 
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APPENDIX B 


Algorithms for the Complete Elliptic Integrals 

There are a number of existing algorithms which are capable of 
providing extremely accurate results for. the three types of elliptic 
integrals. Because of the frequent number of times it is necessary 
to compute these integrals it was also important to obtain methods 
which were relatively quick. It was for this reason that the below: 
algorithms were chosen for the first and second elliptic integrals over 
methods \diich employed Landen or Gaub transformations*. 


A. Algorithm for the Computation of the First 

-33 

Complete Elliptic Integral 


K(k) = (((0.032024666 k^^ + 0.054544409) k_^ + 0.097932891) k ^ •*- 
1.3862944) - (((0.010944912 k ^ + 0.060118519) ^ 

0.12475074) + 0.5) In (k^^) (Bl) 

vdiere k^^ = 1 - k^ 


B. Algorithm for the Computation oF the 
Second Complete Elliptic Integral^^^ 

E(k) = (((0.040905094 kj + 0.085099193) k^^ + 0.44479204) k^^ + 1.0) 

c c c 

(((0.01382999 kj^ + 0.08150224) + 0.24969795) k^^) In (k^^) 

O ^ CO* 

vfliere k ^ = 1 - k^ (B2) 
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C. Algorithm for the Computation of the 

: ' ■ — 35 — 

Third Complete Elliptic Integral 

Existing algoritfins for the coii^tation of the third COT?>lete 
elliptic integral are more conplex than the above algorithms. The method 
of computation en 5 >loyed herein uses the Landen transformation and was 
found to provide excellent results. 
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IF 


U: a = MO X 

n = a/d 

b = 2x(£xn-*'b) 
f = c 
d * n + d 
n =* MD 
c = b/d + c 
MO = + MO 

|n - k^l > CA. X n: 

k^ = 2 
c 

Go to LI 


IF jn - k 1 < CA X nt 

n(Y^, k) = FIFO X c/MO 
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APPENDIX C 


Calculation of the Initial Values for the 
Geometry of the Jet 

For small values o£ the parameters Cj, a, and x, results of 
Spence’s^ two-dimensional theory was utilized for the determination of 
initial values for the jet geometry in the plane of symmetry. The 
necessary equations are summarized below. 
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Because Spence anployed assumptions inherent to linear aerodynamic theory, 
the relation for the curvature, Equation C3, was obtained from the approxi- 
mation 

K i 

dx'^ 
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APPENDIX D 


Curvature Integration Algorithm 

The conqputations of the iteration scheme enqployed for the solution 
of the jet problem are initiated with a set of values for the jet curva- 
ture at the control points. Ihese values are utilized to determine the 
resulting jet slope values at the same locations. This task is not as 
straightforward as might be presmed since the integrals of the curva- 
ture must be properly bounded in order to ensure the evaluation of the 
jet sieves. The definition of curvature 


d 
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along with the condition 


0 ) 1 ) 


(%) x-O ” ■ 

provides a first order differential equation for the jet slope which has 
the solution 
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n 

^ , -tan Ct >c.) _ f pjj 

Jl + tan^ (j * tt) 0 

Obviously and therefore the integral of the cxirvature must be properly 
bounded in order to ensure that the term under the radical of Equation 
(D2) does not become negative (|Aj £ !)• addition, it is expected 
that the slope of the jet remains negative. Unfortunately, the set of 
curvature values predicted from a previous iterate may not satisfy this 
criterion. This dilemma is circumvented by bounding the integral 

CO 00 

I k' dx* = I k' dx' I k' dx' (D4) 

0 0 

where x^ is the largest control point. 

Because the curvature exhibits a logarithmic behavior, it has been 
established that the best method for obtaining the integral of the 
curvature at each of the N control points is by analytically integrating 
the approximation 



The coefficients fb ) in the above series are determined for each interval 

m 

- Xj^) by satisfying Equation (D5) at the M+1 control points closest 
to this interval. In this manner a set of values for the integrals 
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necessary for the evaluation of (cf . Equation D3) and ultimately the 
first integral on the right hand side of Equation (D4) may be determined. 

The remaining integral of Equation (D4) is obtained by anploying 
the series approximation for the slope (refer to Section A2 of Chapter 
III). 
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For large values of x, the above series may be truncated to provide the 
appr oxiinat ions 
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Substitution of Equations (D7) and (D8) into (Dl) provides the resulting 
approximation for the curvature 
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Bnploying Equation (D9) , it is now possible to obtain an approxiiiate 
value for the remaining integral of Equation (D4) 
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A value for the right h^d side of Equation (1)11) has already been obtained. 
Therefore, by assuming a value for (dfj/dx)^ (one would expect the slope 

ratio of the values of the left and right hand sides of Equation (Dll) 
may be obtained. This ratio provides a multiplying factor for the set 
of integral values of Equation (D6) resulting in a properly bounded set. 


of the jet at infinity to lie within the range 
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APPENDIX E 


Calculation of the Grid Point Locations 

The region of the jet located closest to the airfoil obviously 
supplies the dominant effects on the wing not only because of its 
proximity but also because it contains the largest values of the vortex 
strength. It is therefore advantageous to employ more control points 
in this region, the number increasing with increasing values of Cj, x, 
and a. For this rejison, the grid points (vMch ultimately provide the 
control point locations) are chosen to lie at equidistant intervals along 
the jet sheet. 

With the 'jet ;5lope known as a continuous function of x (refer to 
Section A2 of Chapt<jr III) the arclength (K) along the jet is determined 
from the first order, linear ordinary differential equation 


d5 

35 ^ 



(X) 


(El) 


where 


5(x » 0) » 0 


The solution to this* system is obtained by utilizing a fourth order Runge- 
23 

Kutta technique to obtain the necessary starting values and a fourth 
order Adams -Moulton^^ scheftie for the deteimination of the ronaining values. 
For the problem of concern these two techniques may be written as 
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«n*l = ^ [HCV " 2 * VJ 

where A„ = 

and 


^+1 " ^ * 24 “*^^+1^ 


+ 19 B(Xjj) 


5 HCx^.p ♦ E(i^_2J1 


CE2) 


CE3) 


where A = (x - x ) *= constant; 

n+1 n-' 

respectively. 

Given a suitable interval (A? = h), values o£ (i = 1 ,..., 201 
herein) which provide equally spaced values' for the arc length (Cj.) inay 
now be determined. Since y^2 = x^^ = 0) is less than A? the procedure 

is initiated by. choosing X2 = h (the barred variables are used in or 
calculated from Equations (E2) and (E3) , whereas the lonbarred variables 
are those being sought) and calculating the resulting ^2 Equation 
(E2) . If ?2 (refer to Figure 28a) is greater than K2 ~ Ki * (vdiicli 
must be the case for the initial attempt) Xj and are set ociual to X2 
and X2 respectively, t'ne values of X2 and h are halved, and a new value for 
C*2 is determined. This procedure is continued imtil < ^2» is obtained 
(refer to Figure 28b). At this pomt the values of f’2, and T3 •’s well 
as their absissa (x^^, X2, and x^), and ^2» sre known therefore jm iiiter- 
polation scheme may be utilized to obtain X2. 

Having determined the three necessary starting values X2, x^) 
the implicit Adams -Nbul ton technique is now employed.' is first checked 
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A. 


to insure its value is less than 5^ + AS)* If it is not 5^, 

viiich. is the value of T at + h (where h is now the final result of 
the halving procedure described in the previous paragraph), is obtained 
from Equation (J5) . This inarching technique is continued until is 
determined such that The value of is then obtained from the 

interpolation schane (refer to Figure 28c) . This method is continued 
until an array containing the desired number of abscissa (x^) has been 
determined. 

24 

A second order, divided difference interpolation technique was 
chosen for the determination of the abscissa C^Cj^)* This schone has the 
form 

“ ^f^i-2^ * ^V^i-2^ ^i-1^ ^V^i-2^ ^V^i-1^ ^^^i-2» ^i-l» ^i^ 

where x[^. is a divided difference. Therefore, 

X J X 


2 * - 1 ^ ~ 


^ i -1 - ^ i -*2 

H-1 ■ 4-2 




^i - ^i-2 


4 ' ^i-i 4-1 ~ 4-2 
- 5'i-l ’ 


The values I"., C- i» 3jad ^ are chosen such that only is grccUter 

X X“X X*^ X 

than Cjj (refer to Figure 28c). 
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APPENDIX F 


Evaluation of the Integrals Associated With the Jet Problem 
in a Neighborhood of the Control Point 


Consider the problem of evaluating the Cauchy principal value 


x+e- 

f 


x-e. 


x-x' 


dx» 


(FI) 


which involves a singularity at the point x* = x. The function f(x’) 
is continuous and non-singular in the region of concern (refer to 
Section A1 of Chapter IV), Consequently, for small values of and €2* 
a sufficient approximation for this integral may be obtained by en^loying 
a fourth degree polynomial approximation for f(x')» Ihe resulting 
integration may then be performed analytically. 

The method used herein employs Lagrange's interpolation foimila^^ 
for approximating ffx')in the region x^ £ x* < x^ (where x^ < x-Cj and 
X4 > x+9 

4 

£(x') = E Qi(x') f(x.) (F2) 

i=0 ^ ^ 


idiere 



(x' rxp (x' rXj^) (x'. rX^ (x' -Xj^) 


; j ,k,m,n ^ i 


(F3) 


loe 



and 


i 


xi » x-e^ 
X2 « X 

X3 = x+62 

Defining 

Equation CF3) may be written. 


; j ,k,m,n / i (F4) 


Qi - (x' - a^^ X* + a 2 |^ x' - a^^ x* + aQ^ (F5) 

\Aere 

aii = X. Cxj^(x^+x^) + x^ x^) + xj^ x^ 

^2i “ 

for j,k,m,n 7* i. Employing tlie identity x’ « x - (x-x')» Equation (FS) may 
be expressed in the form 

{(x-x')^ •*• [ajj^ - 4x] (x-x')^ + (a 2 j - 3x 6x^) 
Cx-x*)^ + - 2x a 2 ^ + 3x^ a^^^ - 4x^] (x-x*) + 

f^Oi ■ * ®li ®2i " ®3i * 
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Substitution of Equation (F6) into Equation (F2) provides the result 

^ i 

£(x*) = Z B.(x-x‘) 

i=0 ^ 


)diere 


Bq “ f (Xj) A. [3qj - X a^ + X &2 j - X aj^ + x ] 

4 2 3 

B j ® 2 f (Xj ) Aj J ” ^2 j ^ 


3=0 

4 

2 

3=0 

4 


B 2 * f C^jD ^3j ^ 


B, « 2 £(x.) A. [a,. - 4x] 
^ j=0 ^ ^ ^ 


B. - 2 £(x.) A. 
^3=0 ^ ^ 


This last expression for £(x’) may now be used in Equation (FI) to ob 
tain desired approximation 



x-e 


1 



®o 


n 

Cyy 


(F8) 
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TABLE 


COMPUTED RESULTS 


X = 31.3 degrees 



a Cdegrees) 

<x 

<T) 

0 (degrees) 

S 

2.08 

0 

2.548 

-1.777 

8.62 

0.02874 

2.08 

5 

3.259 

-1.499 

11.01 

0.04397 

2.08 

10 

3.925 

-1.113 

13.28 

0.06248 

2.0 

0 

2.434 

-1.709 

8.22 

0.02662 

3.0 

0 

3.326 

-2.563 

10.37 

0.04341 

4.0 

0 

4.197 

-3.418 

12.35 

0.06199 

6.0 

0 

5.817 

-5.127 

15.43 

0.09904 

8.0 

0 

7.388 

-6.836 

18.00 

0.13562 
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a) Jet Flap. 




Figure 1. Powered Flap Configurations, 
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Figure -2. Jet -Flap and Airfoil Geometry, 
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Figure 8. Horseshoe Vortex System. 
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b) Location of the Control Point 
on the i^^ Element. 


Figure 9. Finite Element Geometry. 
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Figure 10. Variation o£ the Lift Coefficient vdth Jet 

Monentun Flux for *a ■ 0*. 




RECTANGULAR WING 
FULL SPAN FLAP 
AR=6.8 
Cjs 2.08 

T * 31.3® 


o WILLIAMS a ALEXANDER (EXPT.) 

- MASKELL a SPENCE (THEORY) 

- TOKUDA (THEORY) 

- LISSAMAN (THEORY) 

- PRESENT THEORY 


0 (DEGREES) 


Figure 11. Variation of the Lift. Coefficient vdth 
Angle of Attack’ for Cj = 2.08. 
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Figure 12. 


Variation of the Drag Coefficient with 
Angle of Attack for C, = 2,08, 
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Figure 13, Variation of the Drag Coefficient 
with Jet Momentum Flux, 
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Fijjure 14. Spanwise Distribution oi 
Across the "Airfoil foa 
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Figure 19. Streanadse Distribution 
Across the Airfoil £ 
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Figure 23. Variation o£ the Jet Vortex Strength in the Plane of 
Symmetry vdth Changes in the Angle of Attack. 
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Figure 25. Variation of the Downwash Angle vdth 
Angle of Attack- for C, = 2.08* 



Figure 26, Variation of the Pressure Coefficient with 
Jet Momentum Coefficient for a = 0*. 














